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Abstract. Volcanic SO 2 column amount and injection height 
retrieved from the Ozone Monitoring Instrument (OMI) with 
the Extended Iterative Spectral Fitting (EISF) technique are 
used to initialize a global chemistry transport model (GEOS- 
Chem) to simulate the atmospheric transport and lifecycle of 
volcanic SO 2 and sulfate aerosol from the 2008 Kasatochi 
eruption, and to subsequently estimate the direct shortwave, 
top-of-the-atmosphere radiative forcing of the volcanic sul- 
fate aerosol. Analysis shows that the integrated use of OMI 
SO 2 plume height in GEOS-Chem yields: (a) good agree- 
ment of the temporal evolution of 3-D volcanic sulfate dis- 
tributions between model simulations and satellite obser- 
vations from the Moderate Resolution Imaging Spectrora- 
diometer (MODIS) and Cloud- Aerosol Lidar with Orthog- 
onal Polarisation (CALIOP), and (b) an e-folding time for 
volcanic SO 2 that is consistent with OMI measurements, re- 
flecting SO 2 oxidation in the upper troposphere and strato- 
sphere is reliably represented in the model. However, a con- 
sistent (~25 %) low bias is found in the GEOS-Chem sim- 
ulated SO 2 burden, and is likely due to a high (~20 %) bias 
of cloud liquid water amount (as compared to the MODIS 
cloud product) and the resultant stronger SO 2 oxidation in 
the GEOS meteorological data during the first week after 
eruption when part of SO 2 underwent aqueous-phase oxi- 
dation in clouds. Radiative transfer calculations show that 


the forcing by Kasatochi volcanic sulfate aerosol becomes 
negligible 6 months after the eruption, but its global aver- 
age over the first month is — 1.3 Wm -2 , with the majority 
of the forcing-influenced region located north of 20° N, and 
with daily peak values up to — 2 Wm“ 2 on days 16-17. Sen- 
sitivity experiments show that every 2 km decrease of SO 2 
injection height in the GEOS-Chem simulations will result 
in a ~25 % decrease in volcanic sulfate forcing; similar sen- 
sitivity but opposite sign also holds for a 0.03 pm increase of 
geometric radius of the volcanic aerosol particles. Both sen- 
sitivities highlight the need to characterize the SO 2 plume 
height and aerosol particle size from space. While more re- 
search efforts are warranted, this study is among the first to 
assimilate both satellite-based SO 2 plume height and amount 
into a chemical transport model for an improved simulation 
of volcanic SO 2 and sulfate transport. 


1 Introduction 

Atmospheric sulfate aerosols play an essential role in atmo- 

spheric radiative transfer by scattering and absorbing solar 

and terrestrial radiation (Hansen et al., 1978; Toon, 1982). 
Although tropospheric sulfate aerosols have a lifetime of a 

few weeks, their impact on climate change is persistent and 
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Table 1. Volcanic SO 2 plume modeling studies for the 1991 Pinatubo eruption. 


Studies 

SO 2 plume heights 

Data Sources 

Models 

Boville et al. (1991) 

21-24 km 

Ground-based report 

3-D CCM2 

Tie et al. (1994) 

19-26 km 

SAGE II 

2-D chemical-dynamical-radiative coupled 

Bekki and Pyle (1994) 

21-28 km 

SAGE II 

2-D chemical-dynamical-radiative coupled 

Zhao et al. (1995) 

20-30 km 

SAGE II, Lidar 

1-D chemical 

Timmreck et al. (1999) 

19-27 km 

SAGE II, Lidar 

3-D ECHAM4 

Pitari and Mancini (2002) 

1 8-25 km 

SAGE II 

3-D CCM-CTM coupled 

Liu and Penner (2002) 

19-26 km 

SAGE II, Lidar 

3-D DAO GCM 

Savarino et al. (2003) 

20.5-3 1km 

SAGE II 

2-D chemical-dynamical-radiative coupled 

Lohmann et al. (2003) 

21.5-29 km 

SAGE II 

3-D ECHAM4 


significant (Seinfeld and Pandis, 2006; IPCC, 2007; Wang et 
al., 2008). In contrast, stratospheric sulfate aerosols result- 
ing from volcanic eruptions can have lifetimes of 1-3 yr, and 
hence have more distinct but irregular (or sporadic) effects 
on global atmospheric chemistry and Earth’s radiative en- 
ergy budget (Budyko, 1977; Hofmann and Solomon, 1989; 
Deshler et al., 2006). Indeed, oxidation of volcanic SO 2 gas 
by OH and H 2 O 2 is a major pathway for producing strato- 
spheric sulfate aerosols, which are highly scattering and in- 
crease planetary albedo in the UV and visible, and this in 
turn leads to radiative cooling of the Earth’s troposphere 
and surface (Robock, 2000). Through heterogeneous reac- 
tions, volcanic sulfate aerosols may also affect chlorine (such 
as CIO) and nitrogen (such as HNO 3 ) chemical cycles in 
the stratosphere, impacting ozone production and destruc- 
tion mechanisms (Hofmann and Solomon, 1989; Russell et 
al., 1996; Solomon, 1999). Solomon et al. (2011) showed 
that stratospheric aerosols have increased in abundance in the 
last decade, likely due to a series of moderate volcanic erup- 
tions (Vernier et al., 2011), resulting in a radiative forcing of 
~ —0.1 Wm “ 2 in average, counteracting the positive forcing 
due to anthropogenic CO 2 . 

Since the Nimbus 7 Total Ozone Mapping Spectrom- 
eter (TOMS) detected the SO 2 clouds from the El Chi- 
chon eruption in 1982 (Krueger, 1983), satellite measure- 
ments have been an indispensable tool for characterizing the 
spatio-temporal distribution of global volcanic SO 2 emis- 
sions. These measurements use the strong SO 2 absorption 
band at 305-330 nm for retrieval of the total column amount 
of SO 2 . Cam et al. (2003) derived a long-term record of vol- 
canic SO 2 emissions from the TOMS satellites that were 
in operation nearly continuously from 1978 to 2005. This 
data record is now being continued and improved with the 
OMI data (Levelt et al., 2006), and supplemented by so- 
lar backscatter ultraviolet (SBUV) measurements from the 
Global Ozone Monitoring Experiment (GOME) (Burrows 
et al., 1999) and GOME-2 (Munro et al., 2006), and by 
infrared (IR) measurements from the Moderate Resolution 
Imaging Spectroradiometer (MODIS) (Watson et al., 2004), 
Atmospheric Infrared Sounder (AIRS) (Prata and Bernardo, 
2007), Infrared Atmospheric Sounding Interferometer (IASI) 


(Karagulian et al., 2010) and Advanced Spacebome Thermal 
Emission Spectrometer (ASTER) (Pieri and Abrams, 2010). 

However, while satellite-based SO 2 emission inventories 
provide climate models with a unique description of the 
spatio-temporal distribution of volcanic SO 2 , they provide 
limited information on the SO 2 vertical distribution. Con- 
sequently, current practice is to specify the SO 2 injection 
height using the volcanic explosivity index (VEI) of the erup- 
tion, which is assigned based on many observable parameters 
available from ground-based reports and is not necessarily an 
accurate indicator of volcanic SO 2 injection height (Spiro et 
al., 1992; Simkin and Siebert, 1994; Andres and Kasgnoc, 
1998; Robock, 2000). Indeed, in the TOMS SO 2 retrieval al- 
gorithm, the SO 2 is assumed to be homogenously distributed 
below either 5 km or 20 km altitude (Krueger et al., 1995, 
2000). Hence, the lack of observation-based characterization 
of SO 2 plume height has led to various discrepancies in quan- 
tification of the climatic effect of volcanic sulfate aerosols. 
As an example, Table 1 shows a list of different volcanic SO 2 
plume heights used in various modeling studies of the 1991 
Pinatubo eruption; some studies derived the SO 2 injection 
height based on the same SAGE (The Stratospheric Aerosol 
and Gas Experiment) aerosol product, but obtain different es- 
timates. We note that volcanic sulfate aerosols are the result 
of oxidation of SO 2 , and hence, the difference in gravita- 
tional settling velocity of SO 2 gas and aerosol particles as 
well as the vertical variation of atmospheric oxidation ca- 
pacity can yield discrepancies between the shapes of vertical 
profiles of volcanic aerosols and SO 2 plumes. Hence, while 
aerosol vertical profiles can be a good proxy for SO 2 plume 
injection heights, the direct retrieval of SO 2 plume height 
from satellite measurements is highly advantageous and is 
expected to improve modeling of the temporal variation and 
climatic effects of volcanic aerosols. 

In recent years, new-generation satellite instruments such 
as the polar-orbiting hyperspectral UV sensors (e.g., OMI) 
and advances in retrieval techniques have expanded our abil- 
ity to measure volcanic emissions (Clarisse et al., 2008; Eck- 
hardt et al., 2008; Yang et al., 2009, 2010; Rix et al., 2012) 
beyond the total SO 2 column amount. In particular, Yang 
et al. (2009, 2010) developed an extended iterative spectral 
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fitting (EISF) technique to simultaneously retrieve both SO 2 
amount and SO 2 altitude from OMI measurements. They 
found that EISF retrievals of SO 2 plume height were in good 
agreement with other observations, and their estimate of SO 2 
amount has higher accuracy than those derived from the (op- 
erational) OMI linear fit retrieval algorithm (Yang et al., 
2007). 

To demonstrate the value of these advances in remote sens- 
ing of SO 2 plumes for climate studies, in this paper we 
use EISF SO 2 column and altitude retrievals (as in Yang et 
al., 2009, 2010) to constrain a 3-D global chemical trans- 
port model (CTM; GEOS-Chem) simulation of the volcanic 
aerosol distribution and direct radiative forcing following 
the August 2008 eruption of Kasatochi (Aleutian Islands). 
Kravitz et al. (2010) illustrated the importance of 2008 
Kasatochi volcanic aerosol forcing on a regional scale, al- 
though the climate effect on a global scale appeared insignif- 
icant; they assumed a total SO 2 emission of 1.5 Tg which 
was evenly distributed in three model layers (10- 16 km) of 
a GCM. This study differs from prior modeling studies in 
that: (a) the CTM is initialized with the direct retrieval of the 
amount and injection altitude of volcanic SO 2 from OMI, 
(b) the CTM results are evaluated, and likely causes of un- 
certainties in the simulation of the volcanic SO 2 lifecycle 
from transport to sink terms in the atmosphere are diagnosed, 
with data from multiple A-Train satellite sensors including 
MODIS aerosol products, MODIS cloud products, additional 
OMI SO 2 data that are not used to initialize the CTM simula- 
tion, and aerosol extinction profiles from the Cloud- Aerosol 
Lidar with Orthogonal Polarization (CALIOP), (c) a sensitiv- 
ity study is conducted to analyze the volcanic aerosol forcing 
as a function of SO 2 injection height specified in the CTM. 
We describe the satellite data in Sect. 2, the configuration of 
the GEOS-Chem CTM and the method for calculating vol- 
canic aerosol radiative forcing in Sect. 3, present results of 
the baseline simulation in Sect. 4 and sensitivity simulations 
in Sect. 5, and finally summarize the paper in Sect. 6. 


2 Satellite data 

SO 2 data retrieved from OMI with the EISF algorithm (Yang 
et al., 2009, 2010) are used in this study to initialize and 
validate the SO 2 distribution in the model. The EISF tech- 
nique takes full advantage of the hyper-spectral BUV mea- 
surements from OMI to improve the accuracy of SO 2 column 
retrievals and simultaneously determine the effective altitude 
of the SO 2 plume. It was designed to address the following 
two disadvantages of earlier algorithms: (a) a priori assump- 
tion of the SO 2 vertical distribution that sometimes results in 
large errors in the retrieved SO 2 amount; (b) underestimation 
of SO 2 burdens, especially during large eruptions, because 
the relationship between BUV radiance and atmospheric SO 2 
column increments is assumed to be linear whereas it actu- 
ally becomes non-linear as the SO 2 burden increases. In the 


EISF algorithm, the SO 2 vertical distribution is assumed to 
be Quasi-Gaussian (with a fixed half width of 2 km in this 
study), and hence, the EISF retrievals provide an SO 2 amount 
and effective plume altitude for each OMI footprint (24 km 
x 13 km at nadir). 

The MODIS aerosol optical depth (AOD) product is used 
to validate our simulation of volcanic sulfate aerosol. The 
MODIS instruments aboard NASA’s Terra and Aqua satel- 
lites provide near daily global coverage at their local equato- 
rial overpass times of 10:30 a.m. and 1:30p.m., respectively 
(Remer et al., 2005). Since MODIS AOD is a columnar quan- 
tity that has limited information about the aerosol chemical 
composition and aerosol vertical distribution, a direct com- 
parison between MODIS AOD and the modeled volcanic 
sulfate AOD is not straightforward, in particular when other 
types of aerosols dominate in the atmospheric column. How- 
ever, over remote regions where background AOD is gen- 
erally low, the spatial distribution of high MODIS AOD is 
still expected to be a good indicator of the transport path or 
distribution of volcanic aerosol. Hence, we use MODIS AOD 
for the evaluation of model- simulated transport pathways and 
distributions (instead of the absolute amount) of volcanic sul- 
fate aerosol. For this purpose, we use the MODIS level 3 
AOD product (from both Terra and Aqua) with a spatial res- 
olution of 1° x 1° and an uncertainty of ±0.05 AOD ±0.03 
over the ocean and ±0.20 AOD ±0.05 over the land (Remer 
et al., 2005). 

Since in-cloud oxidation is a major sink for volcanic SO 2 , 
the MODIS(MOD08) level 3 cloud product (King et al., 
2003) is used to evaluate the accuracy of cloud liquid wa- 
ter and cloud fraction in the GEOS-Chem model, which then 
provides a basis for the interpretation of any differences be- 
tween the GEOS-Chem simulated and OMI-observed SO 2 
distribution. The cloud information in the MOD08 product is 
a result of 1° x 1° averaging of MODIS level 2 (MOD06) 
cloud products that include information on cloud particle 
phase (water, ice, or mixed), cloud fraction, and cloud op- 
tical properties (optical thickness and size for both water 
droplets and ice crystals). MOD06 is retrieved through a 
series of algorithms including the cloud mask algorithm at 
1 km resolution (Ackerman et al., 1998; Frey et al., 2008) 
and cloud microphysical retrieval algorithm at 1 km resolu- 
tion (day-time only) (Platnick and King, 2003). This study 
uses the cloud fraction and cloud liquid water path (LWP) 
saved in the MOD08 daily product to compare with counter- 
parts in the Goddard Earth Observing System (GEOS) mete- 
orological fields used by GEOS-Chem. Past analysis showed 
that the biases in the liquid water path data retrieved from 
MODIS depend on the cloud type, 3-dimentinoal structure 
of the clouds (e.g., broken vs. overcast, adiabatic vs. non- 
adiabatic, etc.), satellite-Sun-Earth geometry, and whether or 
not having drizzle in the clouds or absorbing aerosols above 
the clouds (Wilcox et al., 2009; Seethala and Horvath, 2010; 
Min et al., 2012). In the middle-to-high latitude oceanic re- 
gion of our interest, Seethala and Horvath (2010) found that 
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the MODIS liquid water path data overall overestimates the 
counterpart retrieved from space-borne microwave (AMSR- 
E) instrument, although significant underestimation can also 
occur especially for broken clouds. While quantifying the un- 
certainties in MODIS liquid water path product in our study 
region and time period is challenging, all past studies sup- 
port that summation or averaging of MODIS liquid water 
path over a large spatial domain often reduce the uncertainty 
(Seethala and Horvath, 2010; Min et al., 2012), which is also 
the strategy used in this study during the intercomparison of 
MODIS and GEOS-5 liquid water path (Sect. 4.1). 

To evaluate the model simulation of volcanic aerosols in 
the vertical direction, we compare model results with data 
from the Cloud- Aero sol Lidar with Orthogonal Polarization 
(CALIOP) instrument, aboard the Cloud- Aerosol Lidar and 
Infrared Pathfinder Satellite Observation (CALIPSO) satel- 
lite launched in 2006. CALIOP is a two- wavelength (532 
and 1064 nm), polarization- sensitive (at 532 nm) lidar that 
measures atmospheric backscatter with a single-shot verti- 
cal and horizontal resolution of 30 m and 333 m, respectively. 
An extinction-to-backscatter ratio, also referred to as lidar ra- 
tio, is needed to convert the aerosol backscatter to extinction. 
The CALIPSO aerosol algorithm selects a “best-match” lidar 
ratio after a series of steps. ( 1 ) a cloud aerosol discrimina- 
tion (CAD, Liu et al. 2009) algorithm based upon probabil- 
ity distribution functions (PDLs) of layer averages of 532 nm 
backscatter, attenuated total color ratio, the midlayer altitude 
z, and the depolarization ratio is used to separate clouds from 
aerosols, and to differentiate layers of non-spherical dust par- 
ticles from layers of spherical particles (e.g., liquid sulfate); 
(2) based upon the geolocation and season of CALIPSO ob- 
servations as well as the CAD in step (1), the aerosol type and 
the lidar ratio are selected from a look-up table that is gener- 
ated from cluster analysis of AERONET data and in situ ob- 
servations (Omar et al., 2005; Winker et al., 2009; Winker et 
al., 2010). To fulfill feature finding and layer classification re- 
quirements, the current CALIOP level-2 version 3 algorithm 
yields an aerosol profile product at a horizontal resolution of 
5 km and vertical resolution of 60 m under 20 km altitude. 
In this study, the quality control flag in the CALIOP level-2 
product is used to ensure high quality CALIOP retrievals of 
aerosol layers for comparing volcanic sulfate aerosols from 
the GEOS-Chem simulations. 


3 Methodology 

3.1 GEOS-Chem model, simulation initialization, and 
sensitivity experiments 

A global 3-D CTM, GEOS-Chem (Bey et al., 2001), is used 
to simulate the evolution of volcanic SO 2 . The model is 
driven by assimilated meteorological data from the GEOS 
at the NASA Global Modeling and Assimilation Office 
(GMAO). In this study, version 9-01-01 (http: //GEOS-Chem. 


org) is used at 2° x 2.5° resolution with GEOS-5 47-level 3- 
hourly meteorological fields (interpolated at every 15 min- 
utes to match the time step in the GEOS-Chem). Convective 
transport in the model is calculated from the convective mass 
fluxes in GEOS-5 meteorological fields (Wu et al., 2007). 
For boundary layer mixing the non-local scheme is used (Lin 
and McElroy, 2010). The wet deposition schemes for water- 
soluble aerosols (Liu et al., 2001) and for gases (Mari et 
al., 2000) are implemented. Dry deposition is based on the 
resistance-in-series scheme (Wesely, 1989), with the consid- 
eration of the hygroscopic growth of aerosol particles (Park 
et al., 2004). Anthropogenic emissions of SO 2 in the model 
use as default the EDGAR 3.2 global inventory for 2000 
(Olivier and Berdowski, 2001). The model also uses global 
biofuel emissions (Yevich and Logan, 2003), anthropogenic 
emissions for black carbon and organic carbon (Bond et al., 
2007), shipping emissions from ICOADS (Lee et al., 2011), 
biomass burning from the GFED-2 inventory (van der Werf 
et al., 2009), and a lightning NO x emissions algorithm (Price 
and Rind, 1992). Eruptive and non-eruptive volcanic SO 2 
emissions for each year are implemented in the model using 
the AEROCOM hindcast emission data (Fisher et al., 2011), 
but for the Kasatochi volcanic emissions, we use the OMI 
EISF data (see description below). The default eruptive vol- 
canic SO 2 data provide daily emissions that are on a generic 
1° x 1° grid and are re-gridded into 2° x 2.5° resolution in 
the model. 

Aerosol simulation in GEOS-Chem includes the sulfate- 
nitrate- ammonium system (Park et al., 2004), carbonaceous 
aerosols (Park et al., 2003), sea-salt (Alexander et al., 2005), 
and mineral dust (Fairlie et al., 2007), and couples with 
gas-phase chemistry (Jacob, 2000) through nitrate and am- 
monium partitioning (Park et al., 2004), sulfur chemistry 
(Chin et al., 1996; Alexander et al., 2009), secondary or- 
ganic aerosol formation (Fu et al., 2008), and uptake of 
acidic gases by sea salt and dust (Evans and Jacob, 2005; 
Fairlie et al., 2010). GEOS-Chem includes all major sink 
terms for SO 2 in the atmosphere, including oxidation by 
the hydroxyl radical (OH) in the gas phase and by ozone 
(O 3 ) and hydrogen peroxide (H 2 O 2 ) in the aqueous phase 
at temperatures above 258 K (Fisher et al., 201 1; Wang et al., 
2008a). Stratospheric chemistry in GEOS-Chem is based on 
climatological representation of species sources and sinks, 
and uses the Linoz algorithm of McLinden et al. (2000) 
to simulate stratospheric O 3 (http://wiki.seas.harvard.edu/ 
geos- chem/index.php/Stratospheric .chemistry). The sulfate 
aerosols are partly or totally neutralized by ammonia (NH 3 ), 
and ammonia and nitric acid are partitioned between the gas 
and the sulfate-nitrate- ammonium aerosol phases using the 
ISORROPIA II thermodynamic equilibrium model (Foun- 
toukis and Nenes, 2007). A good agreement with no sys- 
tematic bias was found at the continental scale for compar- 
ison of the GEOS-Chem simulated distribution of sulfate- 
ammonium particles and their extent of neutralization with 
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those from ground-based observations (Park et al., 2004; 
Martin et al., 2004). 

The current simulation of the evolution of volcanic SO 2 
emitted by the Kasatochi eruption was initialized with the 
spatial distribution of OMI EISF SO 2 amount (Fig. la) and 
effect height (Fig. lb) on 8 August 2008. SO 2 plumes with 
column amounts up to 250 DU and effective altitudes up to 
10km can be seen around 52° N, 165° E (Fig. la and b). 
Blocked by a ridge with center line along 155° W (Fig. lc), 
the plume was unable to move eastwards, but instead circu- 
lated around a low pressure system (centered around 50° N, 
170° W) following the anti-clockwise cyclonic flow, and 
hence quickly diluted westwards, to 50 DU with an effec- 
tive height of 4-6 km in the downwind region around 50° N, 
172° W (Fig. la and b). Based upon the distribution of the 
SO 2 amount and effective altitude respectively in Fig. la and 
b, the vertical distribution of SO 2 (as a function of altitude) 
is computed under the assumption that its shape follows the 
Quasi-Gaussian distribution function with a fixed half- width 
of 2 km (but different effective altitude). This assumption is 
consistent with that in the EISF algorithm (Yang et al., 2010). 
The resultant 3-D distribution of SO 2 mass is re-gridded into 
the GEOS-Chem 3-D grid space to be assimilated into the 
model (Fig. lc). In addition, because OMI only provides a 
snapshot of the distribution of SO 2 during the eruption and 
also likely missed the western-most part of SO 2 clouds in 
the study domain of Fig. 1, the estimate of 1.5 Tg of erup- 
tive SO 2 from OMI retrievals may have a low bias (Yang 
et al., 2010). Consequently, a total of 2.0 Tg SO 2 emission 
is specified with the effective injection height of 10 km at 
the model gridbox for Kasatochi. Based upon Waythomas et 
al. (2010), the eruption duration is assumed to be 24 h (on 8 
August 2008) in the model. 

It is worthy noting that the assimilation of OMI SO 2 into 
the model took place in the hour of OMI overpass time, i.e., 
23:00UTC on 8 August 2008, and the assimilation here es- 
sentially is the replacement of SO 2 field simulated by GEOS- 
Chem with the OMI SO 2 field (e.g., similar as model ini- 
tialization). To avoid the discontinuity of SO 2 field due to 
this replacement in the model, a Barnes smoothing technique 
(Barnes, 1964) is used, in which the influence of the innova- 
tion (e.g., difference between OMI and molded SO 2 ) at the 
model grid box (having OMI SO 2 data) on the change of SO 2 
in another gridbox (not having OMI SO 2 ) is inversely pro- 
portional to the distance between these two grid boxes. The 
purpose of this assimilation is to maximum the use of what 
OMI observed to correct the model simulation that otherwise 
would be fully dependent on the specification of the volcanic 
SO 2 point source function in the model. 

Through chemistry inverse modeling constrained by the 
SO 2 column amount retrieved from AIRS, OMI, and GOME- 
2, Kristiansen et al. (2010) estimated that the Kasatochi SO 2 
emissions may have two peaks at 7 and 12 km above the sea 
level, and some up to 20 km. This estimate, for the bulk, is 
consistent with Fig. 1 that shows the peak of SO 2 mixing 





180° 170°W 160°W 150°W 140°W 


0. 50. 100. 150. 200. 250. (DU) 

Fig. 1. Kasatochi volcanic SO 2 plume effective height (a) and SO 2 
column amount (b) as retrieved by the OMI EISF algorithm at the 
OMI footprint resolution on 8 August 2008 (Yang et al., 2009, 
2010), and the corresponding SO 2 column amount mapped onto the 
GEOS-Chem grid box (c). (c) is used to initialize the SO 2 distri- 
bution in the model. See text for details. The pink solid lines in (c) 
are isopleths (at 25 m intervals) of the 500 hPa geopotential height 
(in m). For illustration purpose, the SO 2 data in (c) is interpolated 
at 1° x 1° resolution, but is gridded into 2° x 2.5° resolution in the 
model simulation. 
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< 0. 20. 40. 60. 80. 100. (DU) < 0. 2. 4. 6. 8. 10. (km) < 0. 20. 40. 60. 80. 100. (DU) -20, -10. 0. 10. 20.(DU) 


Fig. 2. (a) and (c): distribution of SO 2 column amount (in DU) and effective height (EH) (in km) of SO 2 on 9 August 2008 as simulated 
by GEOS-Chem (GComLinit) with model initiation of OMI EISF retrieved SO 2 on 8 August 2012; (b) and (d): the respective counterparts 
of (a) and (c) retrieved from the OMI EISF algorithm; (e): same as (a) but from GEOS-Chem simulations without model initiation of EISF 
retrieved SO 2 (GCno_omi), and (f) shows the difference between (a) and (e). (g)-(l) are respectively the same as (a)-(f) but for 10 August 
2008. (m)-(r) are the same as (a)-(f) but for 11 August 2008. The pink solid lines in (a), (g) and (m) are isopleths (at 25 m intervals) of the 
500 hPa geopotential height (in m). L in red in (a), (g) and (m) shows the location of low pressure systems, while A, B, and C in (g) and (m) 
respectively mark the three different transport pathways for SO 2 . Subscripts omi _init and no_omi respectively denote the simulation with and 
without initialization of OMI SO 2 data. 


ratio is in the range of 6-10 km. However, it also is noted 
that in model simulation by Kristiansen et al. (2010), the 
SO 2 emission is specified for two days at the model grid box 
where Kasatochi is located, and hence their scheme for the 
initialization of emission is similar to this study, although no 
OMI SO 2 retrieval are directly assimilated in their model. 
Nevertheless, both the retrieval of SO 2 amount (such as those 
from standard OMI product) and retrieval of SO 2 height 
(such as from research algorithm developed by Yang et al., 
2009, 2010) have uncertainties with best estimate of 20% 
and 1-2 km respectively; low bias in height retrieval often 
corresponds to high bias in SO 2 amount retrieval, and vice 
versa. To investigate the impact of the SO 2 injection height 
(used after the OMI satellite overpass) on the simulation re- 
sults, sensitivity simulations are conducted with different in- 
jection heights of 2, 4, 6, and 8 km (Sect. 4.2). 

3.2 Radiative forcing calculations and sensitivity 
analysis 

Our forcing calculation follows Wang et al. (2008) but 
with improvement in the treatment of cloud effects. A four- 
stream broadband radiative transfer model (RTM), employ- 
ing monthly-mean surface reflectance data (Koelemeijer et 
al., 2003) and the simulated 3-D aerosol sulfate mass is em- 


ployed for the forcing calculations (Fu and Liou, 1993; Wang 
et al., 2004). The RTM is applied to the solar spectrum for six 
bands, ranging from 0.2 to 4 pm. The GEOS-Chem simulated 
volcanic sulfate mass is converted to AOD following Wang 
et al. (2008) in which the hygroscopic effect on sulfate par- 
ticle size and refractive index is considered. Band averages 
of relative-humidity dependent single scattering properties 
of sulfate aerosols (e.g., single scattering albedo, extinction 
cross section, and asymmetry parameter) are tabulated in the 
RTM for computational expediency, while the cloud optical 
thickness is adopted from the GEOS-5 meteorological field. 
In the RTM calculations, all aerosol and cloud particles are 
assumed to be externally mixed (Wang et al., 2008). The dif- 
ference between upwelling solar irradiances calculated in the 
presence and absence of sulfate aerosols, without (with) con- 
sidering the cloud in the RTM calculation, is the clear-sky 
(all-sky) sulfate direct radiative forcing. In each grid cell, the 
forcing calculation is conducted every 6 hours because the 
input cloud properties have 6-hourly temporal resolution. 

The optical properties of sulfate particles are based on 
Wang et al. (2008) in which the size distribution of sulfate 
particles is assumed to have a lognormal size distribution 
with geometric dry radius of 0.07 pm and standard devia- 
tion of 1.8 pm. While this set of optical parameters is typi- 
cal for tropospheric sulfate aerosols that often occur in the 
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neutralized form of ammonium sulfate (Wang et al., 2008), 
stratospheric volcanic sulfate aerosols may less neutralized 
(more acid) and thus have greater hygroscopicity than tropo- 
spheric sulfate aerosols (Russell et al., 1996). Indeed, within 
3-6 months after Pinatubo eruption, the effective radius (or 
equivalently, geometric mean radius assuming no change in 
geometric standard deviation) of stratospheric aerosols was 
shown an increase by a factor of 2-3 (Russell et al., 1996). 
Wang et al. (2008) estimated that for the same amount of sul- 
fate mass with the same size distribution at RH = 5 %, am- 
monium sulfate, ammonium bisulfate, and sulfate acid par- 
ticles can have 20-30 % difference among the radiative forc- 
ing efficiencies (normalized to sulfate mass) at RH = 80%; 
this difference is primarily due to their different hygroscopic 
growth. To consider the uncertainty due to hygroscopicity 
and other factors (such as particle coagulation that are not 
included in the current GEOS-Chem simulation) in the esti- 
mate of particle size, we conducted sensitivity experiments 
to compute the forcing with different sets of sulfate opti- 
cal properties with increasing particle geometric radius from 
0.07 pm to 0.19 pm (Sect. 3.3). 


4 Results 

4.1 Baseline results for SO 2 and volcanic sulfate AOD 
distribution 

The model simulation shows that the SO 2 plume, as a whole, 
moved toward the east after the eruption (Fig. 2a, g and m). 
The initial SO 2 plume center at 52° N, 165° W on 8 Au- 
gust 2008 dispersed toward the southeast on 9 August 2008 
(Fig. 2a) as a result of the southeastward rotation of the 
major- axis of the low pressure system (originally centered on 
52° N, 168° W on Fig. lc), and on 10 August moved 50° N, 
148° W (Fig. 2e). From this center of SO 2 mass on August 
10 , SO 2 plume extends in several directions (Fig. 2 g): (i) to- 
ward the southwest (e.g., location A in Fig. 2g) as a result of 
the blocking ridge along 130° W, (ii) toward the northeast 
(location C) under the influence of a low-pressure system 
centered around 65° N, 150° W. However, the flow toward 
the northeast bifurcates at 55° N, 135° W, with one branch 
continuing northeast, while another branch moves southeast 
(location B in Fig. 2g) and then turns northwest in the west- 
erly anti-clockwise flow circulating around another low pres- 
sure system centered around 52° N, 125° W. The general 
SO 2 distribution on 10 August (Fig. 2 g) is maintained on 
1 1 August, except that the SO 2 cloud is translated eastward 
by ~10° latitude with its center at 55° N, 135° W, and SO 2 
amounts are more diluted to 20-40 DU on average (Fig. 2m). 

Generally good agreement can be found between the mod- 
eled SO 2 spatial distribution (Fig. 2a, g and m) and the OMI 
retrieved SO 2 amount (Fig. 2 b, h and n), especially in terms 
of the location of the volcanic cloud core on 9 August 2008 
and the bifurcation of the SO 2 plume on both 11 and 12 Au- 


gust 2008. However, the modeled patterns overall are more 
diffuse than the OMI observations (Fig. 2), likely reflect- 
ing the difference between the GEOS-Chem model grid size 
(2° x 2.5°) and the OMI footprint size (24 km x 13 km at 
nadir) and the non-ideality in the model (as discussed below). 
In addition, the inability of OMI to retrieve SO 2 located be- 
neath clouds also can partially explain why the flow of SO 2 
is not as continuous and smooth as that in the model simula- 
tions. 

In addition to the overall agreement between modeled and 
OMI retrieved spatial distribution of SO 2 amount, the GEOS- 
Chem modeled distribution of SO 2 effective height on 9-1 1 
August 2008 (Fig. 2c, i, and o) is also consistent with the 
counterparts of OMI EISF retrievals (Fig. 2d, j, and p). Both 
model and OMI retrievals show that the core of SO 2 plume 
was generally maintained at the effective height of 10 km, 
but the effective height for the part of SO 2 plume in the 
southwest direction (C in Fig. 2g and C in top-right corner 
of Fig. 2m) decent to about 2-4 km. In the comparison, it is 
noted that OMI is not sensitive to SO 2 plumes at low altitude. 

Contrast between GEOS-Chem simulations with and with- 
out initialization of using OMI EISF SO 2 data shows sig- 
nificant differences (up to zb 20 DU) in the spatial pattern of 
SO 2 distribution (Fig. 2f, 1, and r) during 9-11 August 2008. 
However, in comparison to the OMI retrievals, it is clear that 
the simulation with the assimilation of OMI EISF SO 2 data 
gives a better description of the SO 2 transport. The simula- 
tion without assimilation of OMI EISF SO 2 appears to give a 
faster dilution of SO 2 from the core, and therefore, an over- 
estimation of SO 2 in locations such as A, B, and C marked on 
Fig. 2g and 2m (Fig. 2f, 1, and r). In addition, in the southwest 
direction (along location A), the simulated SO 2 without as- 
similation reaches too far to the south (around 30° N) on 10- 
1 1 August 2008, while both OMI and simulation with assimi- 
lation show the SO 2 plume only reaches around 38° N. Quan- 
titatively, the simulated SO 2 with assimilation have a corre- 
lation coefficient of 0.73 and normalized root-mean square 
difference of (with respect to OMI retrievals) 1.25, while 
without assimilation respectively have 0.63 and 1.53, with 
OMI SO 2 (after re-gridding to GEOS-Chem grid). Such im- 
provement through assimilation may be useful to short-term 
prediction of the movement of volcanic SO 2 , which can have 
important implications for aviation weather forecast. 

Quantitative comparisons with OMI SO 2 data show that 
GEOS-Chem underestimates SO 2 columns by over 20 DU 
in the plume core on August 10, and by ~60DU on Au- 
gust 11 (Fig. 2). Further comparison of the 40-day time 
series of modeled and EISF-retrieved total SO 2 burden af- 
ter the eruption shows that the model underestimation ap- 
pears persistent throughout the simulation (Fig. 3). Figure 3a 
also shows that the modeled evolution of total SO 4 - mass 
is consistent with the temporal evolution of total SO 2 mass, 
and is peaked around the end of August and early Septem- 
ber. Heard et al. (2012) showed that with the use of their 
profile 1 of SO 2 emission (i.e., 10 % in 5. 5-6. 5 km, 50% 
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Days after July 31, 2008 



Fig. 3. (a) Time series of daily total volcanic SO 2 mass (in log-scale 
on y-axis) after the Kasatochi eruption. Red diamonds and blue tri- 
angles are the results from the GEOS-Chem simulation and OMI 
retrievals, respectively. Also overlaid is the GEOS-Chem simulated 
total volcanic S0 4 mass (green squares). The solid black line is 
a linear least- squares-fit between the GEOS-Chem simulated SO 2 
mass (log scale) and the number of days after the eruption, from 
which an e-folding time of 9 days for SO 2 is derived. Also shown 
at the top right is the equation for a linear least- squares-fit between 
GEOS-Chem simulated and OMI retrieved SO 2 as well as their lin- 
ear correlation coefficient ( R ). (b) Scattered diagram of daily to- 
tal volcanic SO 2 mass between GEOS-Chem (with each injection 
height) and OMI retrievals. Note, the time series of OMI SO 2 data 
is obtained from Krotkov et al. (2010), in which the data points in 
the first two days are estimated based upon the extrapolation to ac- 
count for OMI’s sampling bias due to limited spatial and temporal 
coverage. 


7.5-12 km, and 40% 12- 14 km), their modeled averages of 
sulfate AOD over the Northern Hemisphere between 0° N 
and 85° N for Kasatochi eruption peaked in early September 
centered around 7 September, which is about 9 days earlier 
than that is derived based upon OSIRIS AOD (at 750 nm) 
data. Figure 3 in Kravitz et al. (2012) further showed that 
the peak value (0.006) of zonal averages of OSIRIS AOD 
first appeared on 1 September at 55° N and then expanded to 
the northern latitude region in the following 30-40 days. It 


is noted that the statistics from OSIRIS AOD data can be af- 
fected by its limitations in spatial and temporal sampling (be- 
cause OISIRIS is a limb sensor measuring the visible light). 
Furthermore, the value of averages of AOD not only depends 
on the S0 4 mass, but also is subject to how S0 4 mass 
is distributed spatially, the aerosol scattering properties and 
the relative humidity simulated or prescribed in the model. 
Hence, instead of using AOD or S0 4 - to quantitatively eval- 
uate the model, we conduct the quantitative comparison with 
SO 2 data derived from OMI, and further qualitatively evalu- 
ate the model with AOD data from MODIS and CALIOR 

Based on the OMI SO 2 data for 14-31 August 2008 (i.e., 
blue triangles in Fig. 3a), Krotkov et al. (2010) estimated 
an e-folding time for the Kasatochi SO 2 of around 9 days. 
Interestingly, an identical e-folding time is obtained in our 
GEOS-Chem model simulation for the same time period, 
which suggests that the oxidation rate for converting SO 2 
to sulfate (e.g., the first-order sink rate) in the upper tropo- 
sphere and low-level stratosphere has no systematic bias in 
the model, and is consistent with that derived from OMI re- 
trievals. However, the persistent underestimate in the mod- 
eled SO 2 columns as shown in Figs. 2 and 3a may reflect 
a larger sink term or overestimation of oxidant abundance 
during early plume evolution when SO 2 underwent aqueous- 
phase oxidation in clouds. This hypothesis is evaluated be- 
low by comparing the cloud fraction and LWP in the GEOS- 
5 meteorological fields with those retrieved by MODIS be- 
cause in-cloud oxidation is a major sink for atmospheric SO 2 . 
Figure 3b also illustrates the comparison of daily total SO 2 
mass between GEOS-Chem and OMI-EISF retrievals. As in- 
jection height is lowered from 10 through 2 km, the reduction 
rate of total SO 2 mass drastically increases with associated e- 
folding time decreasing from 9 days to ~3 days, indicating 
that the life time of volcanic SO 2 more depends on injection 
height rather than injected mass. 

Figure 4a illustrates a cloud fraction comparison be- 
tween GEOS-Chem and MODIS for each cloud fraction bin 
(i.e., 0.1) over the region of SO 2 cloud transport (30-70° N 
and 100-175° W) during the early period after the erup- 
tion (about 7 days). The GEOS-Chem or GEOS-5 meteorol- 
ogy fields have a low bias when compared to MODIS for 
cloud fractions less than 0.7, and high bias when cloud frac- 
tions are over 0.8. In addition, the probability density func- 
tion (PDF) of cloud fraction shows that GEOS-5 fields have 
large (low) frequency of small (larger) clouds when com- 
pared to MODIS (Fig. 4a). The comparison of LWP, a first- 
order indicator of water cloud geometric thickness, between 
GEOS-Chem and MODIS (Fig. 4b) shows a similar pattern 
to that of cloud fraction, indicating that: (a) GEOS-5 under- 
estimates the LWP for thin water clouds and overestimates 
the LWP for thick clouds, and (b) GEOS-5 simulates more 
smaller and thin clouds and fewer large and thick clouds. 
Consequently, the GEOS daily total liquid water amount over 
the study region in the first 7 days after the eruption is al- 
ways (10-20%) greater than that characterized by MODIS 


Atmos. Chem. Phys., 13, 1895-1912, 2013 


www.atmos-chem-phys.net/13/1895/2013/ 


J. Wang et al.: Modeling of 2008 Kasatochi volcanic sulfate direct radiative forcing 


1903 




Fig. 4. Comparison (red line) of Cloud fraction (a) and LWP (b) de- 
rived from GEOS-Chem and MODIS during 7 days after the erup- 
tion over the region (30-70° N and 100-175° W). (c) comparison 
of daily total amount of liquid water derived from GEOS-Chem and 
MODIS for the same time period and region. Also shown in (a) 
and (b) are the probability density function (PDF, right y-axis) for 
each cloud fraction (or LWP) bin from GEOS-Chem (blue line) and 
MODIS (purple line). In each panel, the linear correlation coeffi- 
cient ( R ) and the average difference between GEOS-Chem (GC) 
and MODIS are also provided. 


(Fig. 4c), reflecting the relatively more important contribu- 
tion of (b), i.e., the overestimation of the number of thin 
and small clouds in GEOS-5. Presumably, it is those small 
and thin clouds that can more effectively interact with SO 2 
(because of their large area- to- volume ratio). Therefore, dur- 
ing the early period after the eruption, GEOS-Chem overesti- 
mates the abundance of liquid water clouds (or oxidants) able 
to convert volcanic SO 2 into sulfate aerosols in the simula- 
tion, which partially explains the model underestimation of 
SO 2 as shown in Fig. 3. This is especially likely after further 
consideration that MODIS liquid water path may also has a 
positive bias (that is ~10% in global averages over ocean 
when compared to the AMSR measurements) (Seethala and 
Horvath, 2010). Indeed, our sensitivity experiment shows 
that a reduction of liquid water path by 15 % in the first two 
days in GEOS-5 field results in a 5 % increase in SO 2 to- 
tal amount (7 % and 3 % in 1st and 2nd day respectively). 
It is noted that once SO 2 reaches the upper troposphere and 
stratosphere, its main sink is oxidation by OH, and hence the 
consistency of e-folding time between GEOS-Chem simula- 
tions and OMI observations (Fig. 3) indicates that oxidation 
of SO 2 by OH is well represented in the model. 

For further evaluation of the model, the simulated dis- 
tribution of volcanic sulfate AOD is compared with Terra 
and Aqua MODIS level-2 aerosol data on several selected 
days when the pathway of volcanic sulfate AOD is easily 
discernable (Fig. 5). On August 14 (8 days after the erup- 
tion), MODIS AOD maps indicate that the volcanic sul- 
fate aerosols (with mid- visible AOD >0.5) was mainly lo- 
cated over Alaska, Northern Canada, and Northern Mexico 
(marked respectively as regions A and B in Fig. 5b); such 
distributions are reasonably reproduced by the GEOS-Chem 
model simulation in Fig. 5a. The signature of volcanic sul- 
fate aerosols, two weeks after the volcanic eruption, can 
be further identified over the continental US, Western At- 
lantic Ocean, and Europe (respectively marked as C-E in both 
Fig. 5c and d). A period of about 10 days is enough to trans- 
port the volcanic aerosols not only over the entire Atlantic 
Ocean, but also over Asia. About 17 days after the eruption, 
the volcanic sulfate aerosols are transported zonally from Eu- 
rope to Eastern Asia (e.g., region G, H, and I in Fig. 5e) and 
even meridionally to the Southern Pacific Ocean (region F in 
Fig. 5e). Overall, Fig. 5 shows that the hemispheric distri- 
butions of the volcanic sulfate AOD from the simulation are 
comparable with the MODIS AOD observations. 

In addition to using satellite data for evaluating model- 
simulated column amounts (such as total SO 2 burden and 
AOD), we also use the CALIOP lidar data to evaluate the 
vertical profile of aerosol extinction coefficients simulated 
by the model (Fig. 6). Two CALIPSO orbits (i.e., the solid 
blue lines in Fig. 5a and e) over North America (region B 
in Fig. 5a) and East Asia (region I in Fig. 5e) are used for 
comparison with model simulations. On 14 August 2008, 
CALIOP data indicate a layer (marked as LI and L2 in 
Fig. 6a) with high loading of particles above the tropopause 
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Fig. 5. Spatial distribution of volcanic sulfate aerosol optical depth (AOD) at 550 nm as simulated by GEOS-Chem (left column) and as 
retrieved by MODIS (right column) on 14 (top row), 19 (middle row), and 25 (bottom row) August 2008, respectively. The solid black lines 
in (a) and (e) respectively indicate the ground tracks of CALIOP data that are shown in Fig. 6a and c. The pink rectangles and letters highlight 
the regions of large volcanic AOD, as discussed in Sect. 4.1. 


(black line in Fig. 6b) over North America. Within this layer, 
the CALIOP measurements of depolarization ratio at 530 nm 
and backscattering attenuation at 1062 nm both show nearly 
zero values, and CALIOP layer classification algorithm in- 
dicate that this layer are dominated by aerosols with small 
fraction of cirrus (figures now shown). The model simula- 
tion is able to capture a similar vertical distribution of the 
volcanic sulfate aerosol extinction (marked as LI and L2 in 
Fig. 6b), although the coarser model resolution in the strato- 
sphere cannot resolve the very thin-layer structure of the sul- 
fate aerosols detected by CALIOP data. The stratospheric 
sulfate aerosols are also distinctly detected by CALIPSO on 
the ~1 7 th day after the eruption over East Asia (marked as 
L5, L6, and L8 in Fig. 6c), and those are plausibly captured 
by the GEOS-Chem simulation in Fig. 6d, especially the de- 
scending path of aerosols from Siberia, to East Asia, and to 


the western North Pacific (marked respectively as L5, L6, 
and L8 in both Fig. 6c and d). In addition, CALIOP images 
in Fig. 6a and b also indicate the likely deposition of volcanic 
sulfate aerosols in the middle-to-lower atmosphere such as 
over the south central US (marked as L3 and L4 in Fig. 6a 
and b) and over northeast China (marked as G in Fig. 6c 
and d); these “touch down” features are also seen in the 
similar curtain plots showing the difference in GEOS-Chem 
simulation with and without considering volcanic aerosols, 
although non-volcanic aerosols from local source also con- 
tribute to the high loading of particles in regions L3 and L7 
(figures now shown here). 

Quantitatively, the temporal evolution of daily and zon- 
ally averaged mean sulfate AOD from the model and the 
CALIPSO level 2 aerosol products are illustrated in Fig. 7. 
GEOS-Chem simulations show large amounts of volcanic 
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Fig. 6. (a) vertical distribution of 532 nm total attenuated backscatter (km -1 sr -1 ) measured by the CALIOP lidar, and (b) the corresponding 
distribution of tropopause (black line) and the simulated sulfate aerosol extinction coefficient (km -1 ) at 550 nm for the CALIPSO ground 
track in Fig. 5a on 14 August 2008. (c) and (d) are respectively the same as (a) and (b) for the CALIPSO ground track solid black line in 
Fig. 5e. The pink arrows and letters highlight the regions of large volcanic AOD, as discussed in Sect. 4.1. 


sulfate aerosol in the entire atmosphere (Fig. 7a) with a sig- 
nificant fraction of the aerosol in the stratosphere (Fig. 7b). 
The GEOS-Chem simulated sulfate AOD is compared with 
AOD from CALIPSO in Fig. 7c and d only for altitudes 
over 10 km to minimize the influence of non- volcanic sul- 
fate aerosols in the CALIOP data. It is apparent that the 
model produces a comparable AOD (mean AOD is ~0.06) 


to that observed by CALIOP from 60-80° N. During the first 
30 days after the eruption; the modeled and observed strato- 
spheric AODs are very comparable. The modeled AOD, 
however, is drastically decreased, while the satellite mea- 
sured AOD is only slightly decreased after the first 30 
days. Some discrepancy may be explained by the misclas- 
sification of high cirrus cloud as aerosol in the CALIOP 
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Fig. 7. (a) Temporal evolution of zonally averaged volcanic sulfate 
aerosol optical depth (SAOD) from GEOS-Chem simulations, (b) 
same as (a) but for SAOD above 10 km from the surface, (c) the 
same as (b) but SAOD sampled along the CALIOP ground track, 
and (d) same as (c) but based upon the analysis of CALIOP level-2 
aerosol layer product with the consideration of data quality flags. 


algorithm; indeed, the temporal evolution of the zonally aver- 
aged backscattering ratio from CALIOP as shown in Vernier 
et al. (2011) is more similar to the GEOS-Chem simulation 
in Fig. 7b. In addition, as discussed in Heard et al. (2012), 
volcanic sources other than Kasatochi, such as the ongoing 
eruption from Kilauea volcano (19.4° N, 155.3° W) in 2008 


may also contribute to the slower decrease of AOD above 
10 km in the CALIOP data. About 70 days after the erup- 
tion, stratospheric aerosol is no longer detected in satellite 
data over middle and high latitudes, and the simulated sulfate 
AOD is also very small (~4 x 10 4 ). About 100 days after 
the eruption, the stratospheric sulfate AOD in the model is 
negligible with daily values of 10 -4 . 

4.2 Baseline results for volcanic sulfate forcing 

Figure 8 shows the monthly averaged direct shortwave (SW) 
radiative forcing by the volcanic sulfate aerosols at the top 
of the atmosphere (TOA) after the Kasatochi eruption. In 
monthly and global averages, the clear sky and all sky short- 
wave forcing of volcanic sulfate aerosols are both strongest 
(i.e., most negative) in the first month after the eruption, with 
respective values of —2.0 W m -2 and — 1.3 W m -2 (Fig. Sa- 
fi); they steadily become weaker with respective values of 
— 1.0 and — 0.7 Wm -2 in the second month (Fig. 8e-f), and 
—0.04 and —0.03 in the 6th month (Fig. 8i-j). Geograph- 
ically, the volcanic sulfate aerosols are not transported to 
the Southern Hemisphere until the 2nd month after the erup- 
tion, and do not spread over the whole Southern Hemisphere 
until 30 September 2008. Indeed, most of the volcanic sul- 
fate aerosols remained north of 20° N in the first month af- 
ter the eruption (Fig. 8a). The difference in radiative forc- 
ing between clear and all sky depends, to a large extent, on 
the cloud fraction and relative altitude of the sulfate aerosol 
layer and clouds. Since sulfate particles and cloud droplets 
are highly scattering at visible wavelengths where the solar 
spectrum peaks, cloud layers, whether underlying or overly- 
ing the aerosol layer, generally reduce the spectral contrast 
between the bright aerosol layer and darker land/ocean sur- 
face when viewed from the TOA. As such, cloud layers often 
reduce the clear-sky radiative forcing of scattering aerosols 
(Wang et al., 2008). Hence, comparing the third and fourth 
columns in Fig. 8c, g and k show that the difference between 
clear-sky and all-sky radiative forcing is small (large) in re- 
gions where cloud fraction is low (high) such as over the Sa- 
haran desert and the dry regions of the central-to-east Pa- 
cific (cloudy regions include the southern ocean at 40-60° S, 
the tropics, and high latitudes). As expected, over Greenland 
where surface albedo is high (radiatively acting like a cloud 
layer), the clear and cloudy sky forcing is nearly the same. 

A daily time series of global averages of sulfate aerosol 
all-sky radiative forcing at the TOA is shown in Fig. 9. Soon 
after the Kasatochi eruption, the continuous conversion of 
SO 2 into sulfate results in a steady increase of sulfate AOD 
and a correspondingly stronger SW radiative forcing at the 
TOA. The forcing reaches its peak (— 2.1 Wm -2 ) about 16 
days after the eruption, and thereafter decays exponentially 
with time (Fig. 9). In contrast, Earth Radiation Budget Ex- 
periment data shows that the peak forcing (global average) 
of sulfate aerosols from Pinatubo eruption (15° N, 121° E) 
on 15 June 1991 is at least — 4 Wm -2 in September-October 
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Fig. 8. Averages of volcanic sulfate shortwave (SW) radiative forcing (Wm -2 ) at the top of the atmosphere (TOA) for all sky (a) and clear 
sky (b) conditions, as well as their differences (c) in the first month (a-c), second month (e-g), and the 6th month (i-k) after the Kasatochi 
eruption on 8 August 2008. Also shown are the corresponding distributions of cloud fraction (d, h and 1). Denoted on the top right of each 
panel is the global average (weighted by gridbox area) of the corresponding quantity. 



Days after the eruption 


Fig. 9. Temporal evolutions of the daily, global mean volcanic sul- 
fate radiative forcing (Wm -2 ) at the top of the atmosphere (TOA) 
for 5 different volcanic SO 2 injection heights ranging from 10 km 
(shown in red line) to 2 km (shown in purple) with interval of 2 km 
(see legend on the top right for details). Also shown in the legend, 
corresponding to each injection height used in the GEOS-Chem 
simulation, are days (after the eruption on 8 August 2008) when 
the peak of the shortwave forcing occurs. 


1991 (Minnis et al., 1993). The significant larger effect of 
Pinatubo eruption on climate is in part due to the following 
factors: (a) it ejected ~30Tg of SO 2 up to 30 km above the 
sea level, most of which concentrated in 20-27 km altitude 
(McCormick et al., 1995); (b) at this altitude range over the 
subtropics, the intensity of planetary wave activity and the 
phase of the quasi-biennial oscillation regulates the poleward 
transport, but were shown to be not effective in the North- 
ern Hemisphere in June-July 1992; (c) the spread of SO 2 to 
the subtropics in Southern Hemisphere is found to be unex- 
pectedly faster in June-July 1991, which is attributed to the 
abnormality of the planetary wave activates over the equa- 
tor and southern subtropics (Trepte et al., 1993); (d) conse- 
quently, SO 2 amount was mainly located in the 30° N-30° S 
in the first several months after Kasatochi eruption, and the 
larger and longer solar illumination in tropic and subtropics 
enhance the shortwave forcing of volcanic sulfate particles. 
In comparison, while the volcanic sulfate aerosols from the 
Kasatochi eruption (52.1° N, 175.5° W) can reach the trop- 
ics within 1-2 month, most of them are concentrated in the 
mid- and high-latitudes (due to the westerly waves) during 
Northern Hemisphere fall-winter season (Fig. 8), and their 
amount (1-2 Tg) and e-folding time (9 days) are both much 
less than counterparts (30 Tg and ~1 yr) of Pinatubo eruption 
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(McCormick et al., 1995). Consequently, Kasatochi eruption 
may have only affected the global radiative energy budget for 
about 100 days, and has much less impact on climate. 

The timeline of our simulated volcanic sulfate forcing is 
consistent with model experiment’s by Kravitz et al. (2012) 
that showed all volcanic sulfate aerosols might have been de- 
posited out of the atmosphere by February 2009, and the no- 
ticeable forcing may deceased even quicker. A direct quanti- 
tative comparison with the results from Kravitz et al. (2012) 
is complicated by: (a) the model difference in spatial resolu- 
tion and temporal resolution (4° x 5° and focus of monthly 
scale in their climate model), chemical mechanisms (only 
prescribed OH field available in their model to oxidize SO 2 
in the stratosphere), and cloud fields (simulated solely with 
a climate model), and (b) the definition of sulfate forcing 
in which aerosol forcing feedback on stratospheric thermal 
adjustment are considered in their model. Nevertheless, be- 
cause of sulfate particle is highly scattering (with the single 
scattering albedo value close to 1) in visible and other short- 
wave spectrum, our estimate of global forcing at the TOA and 
surface is very similar, with a global averages of — 1 .3 Wm -2 
in August and —0.7 Wm -2 forcing in September, which ap- 
pear consistent with results in Kravitz et al. (2012) showing 
a — 2 Wm -2 of zonal averages of forcing at the surface over 
the Northern Hemisphere in August and September. 

4.3 Sensitivity experiment to SO 2 injection height 

In order to investigate the effect of volcanic SO 2 plume in- 
jection height on SW radiative forcing, we conducted sensi- 
tivity simulations with injection heights of 2, 4, 6 and 8 km. 
Figure 9 shows that the radiative forcing has a strong depen- 
dence on the volcanic SO 2 injection height. As the injection 
height decreases, the magnitude and the duration of the forc- 
ing decreases. For example, for a 2 km injection height, the 
peak forcing is — 0.6 Wm -2 and the sulfate aerosols influ- 
ence SW radiation for about 35 days, which contrasts with 
— 2.1 W m -2 and about 100 days when injection height is set 
at 10 km. 

Since the temporal evolution of the volcanic sulfate SW 
radiative forcing shown in Fig. 9 appears to follow a lognor- 
mal distribution, the following function is found to provide 
a good fit to the forcing-time curves in Fig. 9 for different 
injection heights: 


Sie -2(lnxS 2 ) 2 



where x = injection height (km), y = parameterized ra- 
diative forcing (Wm -2 ), and the two scale factors are 
a function of injection height x: Si=0.4x + 0.5, and 
S 2 = — 0.00375x + 0.0875. It is found that the parameterized 
SW radiative forcing agrees with the GEOS-Chem simu- 
lated forcing for each injection height as a function of num- 
ber of days after the eruption, with linear correlation coeffi- 


cients generally larger than 0.98 (Fig. not shown). As shown 
in Fig. 10, the parameterization adequately reproduces the 
peak value of the radiative forcing (Fig. 10b) and the timing 
for this peak value (Fig. 10a) for different injection heights. 
Based upon this parameterization, a difference of 2 km in in- 
jection height can lead to a 0.4 % difference in the overall es- 
timate of the forcing effect (e.g., forcing multiplied by time) 
in the whole globe. 

4.4 Sensitivity experiment to aerosol size distribution 

In our baseline simulation, the sulfate particles are assumed 
to have a lognormal size distribution with a geometric radius 
of 0.07 pm and standard deviation of 1.8 pm. In order to in- 
vestigate the impact of volcanic sulfate particle size on SW 
radiative forcing at the TOA, sensitivity experiments were 
conducted to compute the forcing with different sets of sul- 
fate optical properties corresponding to increasing geometric 
radii from 0.07 pm to 0.19 pm with a step size of 0.03 pm. 
Wang et al. (2008) showed that as the particle size increases, 
the associated increase in particle extinction cross section 
outweighs the associated reduction in backscattering, and 
thus results in stronger aerosol forcing. As shown in Fig. 11, 
an increment of 0.03 pm in sulfate particle radius results in an 
enhancement of ~0.l-0.2Wm -2 in the 30-day (after erup- 
tion) and global average sulfate all- sky SW radiative forcing 
at the TOA. 


5 Summary and discussion 

GEOS-Chem, a global chemical transport model, has been 
used in conjunction with constraints from the OMI-EISF re- 
trievals of SO 2 amount and effective height to simulate the 
life cycle of SO 2 and volcanic sulfate aerosols after the 2008 
Kasatochi eruption and to study the resultant impact on di- 
rect shortwave radiative forcing. With the use of the OMI 
EISF-based SO 2 product to initialize the SO 2 distribution in 
GEOS-Chem, the simulated lifetime (with an estimated e- 
folding time of 9 days) as well as the spatial distribution and 
temporal evolution of the volcanic SO 2 burden in the atmo- 
sphere after the eruption are both in good agreement with 
OMI SO 2 observations, suggesting that the oxidation of SO 2 
in the stratosphere (primarily by the hydroxyl radical, OH) 
is reliably represented in GEOS-Chem. However, a consis- 
tent low (~25 %) bias is found in the GEOS-Chem simulated 
SO 2 burden, and comparison with MODIS cloud products in- 
dicates that this is likely due to a high (~20 %) bias in cloud 
liquid water amount and a resultant stronger oxidation of SO 2 
in the GEOS meteorological data during the first week after 
the eruption when part of SO 2 is oxidized by clouds. Further 
evaluation with aerosol products from MODIS and CALIOP 
reveal that GEOS-Chem simulations capture the 3-D trans- 
port pathway of volcanic sulfate aerosols including: (a) their 
longitudinal transport from Alaska to central Canada, and 
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Fig. 10. (a) Scatter plot of the volcanic SO 2 injection height (km) and the days (after the eruption) when the peak of daily and global 
averages of volcanic sulfate shortwave forcing occurs. The data are based upon results in Fig. 9. (b) same as (a) but shows the peak value of 
the daily and global average of volcanic sulfate shortwave forcing. The red lines in (a) and (b) respectively show the results based upon the 
parameterizations as described in Sect. 4.3. 



(/^m) 


Fig. 11. Global and 30-day (after eruption) average of volcanic sul- 
fate shortwave forcing at the top of atmosphere as a function of the 
geometric radius (r g ) used in the log-normal size distribution to de- 
scribe the aerosol optical properties in the forcing calculation. 


then to the south central Great Plains in the first two weeks, 
as well as (b) their zonal transport from high-latitude regions 
of North America to mid- and high-latitude regions in Eu- 
rope and Asia, and the consequent transport from Siberia to 
southeast China within three weeks after the eruption. 

Radiative transfer calculations show that the all- sky direct 
radiative forcing at the TOA due to the Kasatochi volcanic 
sulfate aerosols reached a peak in the late second week and 
early third week post-eruption, with a daily, global average 
value of ~2 W m -2 . Consequently, in global and monthly av- 
erages, the volcanic sulfate forcing from the Kasatochi erup- 
tion peaks up to — 1.3 Wm -2 in the first month after the 
eruption, with majority of the forcing-influenced region lo- 


cated north of 20° N; and then gradually weakens to less 
than -O.lWm -2 four months after the eruption. The vol- 
canic aerosol forcing doesn’t influence the entire Northern 
Hemisphere until the middle of the second month after the 
eruption. It is found that clouds can effectively reduce the 
magnitude of the volcanic sulfate forcing by 20-40 %, on av- 
erage. 

Sensitivity analysis shows that accurate description of the 
SO 2 injection height and the initial 3-D distribution of SO 2 
in the CTM are both critical for reliable simulation of the 
lifetime and spatiotemporal distribution of volcanic SO 2 and 
aerosols after the eruption. For the Kasatochi eruption, it is 
shown that the temporal evolution of the volcanic sulfate 
forcing can be parameterized using a log-normal distribution 
as a function of injection height and number of days after the 
eruption. This parameterization indicates that every 2 km re- 
duction of SO 2 injection height results in a 2 day decrease in 
SO 2 lifetime and 0.4 W m -2 reduction in forcing (in global 
and daily averages). Further sensitivity tests also showed that 
every 0.03 pm increase of geometric particle radius used in 
the log-normal size distribution for describing aerosol opti- 
cal properties leads to ~25 % increase in the magnitude of 
the forcing, although the rate of increase falls off for larger 
geometric radii. 

This study is among the first to assimilate both satellite- 
based SO 2 plume height and column amount into a CTM for 
an improved simulation of volcanic SO 2 transport, which has 
important implications for studies of natural climate forcing 
as well as for forecasts of atmospheric opacity that impacts 
aviation safety. 
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